Astronomy & Astrophysics manuscript no. 7098text 


© ESO 2008 


February 5, 2008 





Letter to the Editor 



Discovery of a flux-related change of the cyclotron line energy in 

Her X-1 

R. Staubert\ N.I. Shakura^, K. Postnov^, J. Wilms^ R.E. Rothschild^ , W. Coburn^, L. Rodina^ D. Klochkov^ 



o 
o 

(N 
0^ 



> 
O 

O 
l> 

o 

O 



' Institut fiir Astronomie und Astrophysik, Universitiit Tubingen, Sand 1, 72076 Tubingen, Germany 
~ Sternberg Astronomical Institute, 119992, Moscow, Russia 

Dr. Remeis Stemwarte, Astronomisches Institut der Universitat Erlangen-Niirnberg, Sternwartstr. 7, 96049 Bamberg, Germany 
* Center for Astrophysics and Space Sciences, University of California at San Diego, La Jolla, CA 92093-0424, USA 
^ Space Sciences Laboratory, University of California, Berkeley, CA 94720-7450, USA 

ABSTRACT 

Aims. We present the results of ten years of repeated measurements of the Cyclotron Resonance Scattering Feature (CRSF) in the 
spectrum of the binary X-ray pulsar Her X-I and report the discovery of a positive correlation of the centroid energy of this absorption 
feature in pulse phase averaged spectra with source luminosity. 

Methods. Our results are based on a uniform analysis of observations by the RXTE satellite from 1996 to 2005, using sufficiently 
long observations of 12 individual 35-day Main-On states of the source. 

Results. The mean centroid energy Ec of the CRSF in pulse phase averaged spectra of Her X-1 during this time is around 40keV, 
with significant variations from one Main-On state to the next. We find that the centroid energy of the CRSF in Her X- 1 changes by ~ 
5% in energy for a factor of 2 in luminosity. The correlation is positive, contrary to what is observed in some high luminosity transient 
pulsars. 

Conclusions. Our finding is the first significant measurement of a positive correlation between Ec and luminosity in any X-ray pulsar. 
We suggest that this behaviour is expected in the case of sub-Eddington accretion and present a calculation of a quantitative estimate, 
which is very consistent with the effect observed in Her X-1. We urge that Her X-I is regularly monitored further and that other X-ray 
pulsars are investigated for a similar behaviour. 

Key words, magnetic fields, neutron stars, - radiation mechanisms, cyclotron scattering features - accretion, accretion columns - 
binaries: eclipsing - stars: Her X-I - X-rays: general - X-rays: stars 



1. Introduction 

■ The X-ray spectrum of the accreting binary pulsar Her X-1 is 
. characterized by a power law continuum with exponential cut-off 
' and an apparent line- like feature, which was discovered in 1975 
; (lTrumperetalJll978l) . This feature is now generally accepted as 
I an absorption feature around 40keV due to resonant scattering 
of photons off electrons on quantized energy levels (Landau lev- 
els) in the Teragauss magnetic field at the polar cap of the neu- 
tron star. The feature is therefore often referred to as a cyclotron 
resonant scattering feature (CRSF). The energy spacing between 
the Landau levels is given by = heB/mgC - 11.6keVBi2, 
where B12 = Z?/10'~G, providing a direct method of measur- 
ing the magnetic field strength at the site of the emission of the 
X-ray spectrum. The observed line energy is subject to gravita- 
tional redshift, z, such that the magnetic field may be estimated 
by B12 = {1 + z)£'obs/ll-6keV. The discovery of the cyclotron 
feature in the spectrum of Her X-1 provided the first ever 'direct 
measurement' of the magnetic field strength of a neutron star, in 
the sense that no other model assumptions are needed. Originally 
considered an exception, cyclotron features are now known to 
be rather common in accreting X-ray pulsars, with more than a 
dozen binary pulsar s now being confirmed cyclotron line sources 
dCoburn et al.ll2002l) . In several obje cts, multiple lines have been 
detected (up to four harmonics, see lStauber3l2003l iHeindl et al.l 
|2004| for reviews). 
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In this Letter we present new results (from the last ten years) 
on the energy of the cyclotron resonance scattering feature in 
the pulse averaged X-ray spectrum of Her X- 1 , we summarize 
our knowledge about its variability, and report the discovery of a 
positive coiTelation between Ec and X-ray luminosity and offer 
a physical explanation. A quantitative description of the physics 
of the accretion column under the condition of sub-Eddington 
accretion is c onsistent with t he obs ervation. Preliminary results 
were given in lStaubert et al.l (l2006l) . 



2. Observations 

Her X-1 is probably the best observed accreting binary X- 
ray pulsar Its X-ray spectrum, including the CRSF, has been 
measured by many instruments since its discovery in 1975 
(iTriimper et al.|[l978h . Here we report on a coherent analysis of 
ten years (1996-2005) of repeated observations b y RXTE, in- 
cludin g a reanalysis of those data already used by iGruber et al.l 
(l200lh . Main-On state (35 day-phases < 0.18) observations 
(where the X-ray flux is the highest) of > lOks duration were 
used. These RXTE measurements are of high quality and yield 
comparatively small uncertainties in the measured quantities. 
The joint spectral analysis of PCA and HEXTE data was per- 
formed using the standard XSPEC/FTOOLS (6.0.4) software. 
For the spectral model we have chosen the highecut model 
which is based on a power law continuum with exponential 
cut-off, and the CRSF is modeled by an absorption line with a 
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Table 1. Recent cyclotron line energy measurements in Her X-1 
by RXTE and INTEGRAL. Uncertainties are at the 68% level. 
The data point from INTEGRAL is f rom Klochk ov et al. ( 20061). 
35 d cycle numbering is according to lStaubert et al.l (Il983h . 



Year 
1990 



Observation 


35 d 


Center 


Line Energy 


max. Flux 


month/year 


cycle 


MJD 


keV 


ASM cts/s 


RXTE 










July 96 


257 


50029.75 


41.12 ±0.55 


7.37 ±0.34 


Sep 97 


269 


50707.06 


40.62 + 0.49 


7.49 ± 0.73 


Dec 00 


304 


51897.69 


40.07 ±0.31 


6.04 ± 0.47 


Jan 01 


305 


51933.67 


39.05 ±0.55 


5.72 ± 0.34 


May 01 


308 


52035.48 


39.93 ± 0.63 


7.15 ±0.50 


June 01 


309 


52071.16 


39.73 ±0.52 


6.93 ± 0.20 


Dec 01 


314 


52245.09 


40.04 ± 0.22 




Aug 02 


321 


52492.96 


40.01 ±0.29 


7.19 ±0.26 


Nov 02 


324 


52599.32 


40.51 ±0.13 


7.64 ± 0.30 


Dec 02 


325 


52634.01 


40.60 ±0.41 


7.55 ±0.34 


Oct 04 


344 


53300.95 


38.51 ±0.51 


4.50 ± 0.24 


July 05 


352 


53577.35 


38.95 ±0.52 


5.12 ±0.37 


INTEGRAL 










July 05 


352 


53576.00 


38.50 ±0.70 


5.12 ±0.37 



Gaussian optical dep th profile . Detai ls of the fitting procedure 
can be found in Cob urn et al. ] (120021) . Here we discuss results 
from spectral analysis of pulse phase averaged spectra only (a 
pulse phase resolved analysis is in progress). Table 1 summa- 
rizes the observation dates, the measured CRSF centroid ener- 
gies and the maximum X-ray fluxes for the corresponding 35 day 
Main-On states, as measured by the RXTE/ASM (formally de- 
termined by fitting the flux history of the respective Main-On 
by a template function derived from the average of many Main- 
On cycles). Table 1 also lists an additional data poi nt from an 
observation with INTEGRAL dKlochkov et al.ll2006l) . 



3. Results 

Fig. 1 displays the measured CRSF centroid energies as a func- 
tion of time, combining the results of our analysis with historical 
measur ements before the RX TE era, as taken from the compila- 
tion by iGruber et al.l (1200 iL their ta bles 2 and 3), as well as the 
2005 measurement by INTEGRAL dKlochkov et al.l2006h . Also 
shown in Fig. 1 (right hand scale) is the maximum flux of the cor- 
responding 35 d Main-On states as measured by RXTE/ASM. 
We report on two main results. 

Firstly, we confirm the apparent difference in the mean cy- 
clotron line e nergy before and after 1991, first pointed out by 
iGruber et al.l (12001.). Taking the measured values of and their 
stated uncertainties at face value, the mean cyclotron line ener- 
gies (Ec) (dotted lines in Fig. 1) from all measurements before 
1991 is 34.9 + 0.3 keV, the corresponding value for all measure- 
ments after 1991 is 40.3 + 0.1 keV (40.2 + 0.1 keV for RXTE re- 
sults only, showing that the very high value measured by BATSE 
is not decisive). However, a comparison of measurements from 
different instruments is difficult because of systematic uncertain- 
ties due to calibration and analysis techniques. Nevertheless, we 
believe that the large difference of ~ 5 keV between the mean 
values and the relative internal consistency within the two groups 
(5 different instruments before 1991 and four after 1991) most 
likely indicate real physics (we will comment on this in the dis- 
cussion). 

Secondly, we present a new result: we have detected a posi- 
tive correlation between the maximum 35 d Main-On state X-ray 



U 



ASM tlux (righl hand m 

Mir HEXE 

HEAO-1 

Balloon HEXE 

GRIS 

Ginga 

BATSE 

BappoSAX 

RXTE 

Integral 



30 I 

42000 



16 ;a 

u 

14 3 



e 
< 



48000 50000 
MJD 



Fig. 1. The centroid energy of the phase averaged cyclotron res- 
onance line feature in Her X-1 since its discovery. Da t a from 
before 1997 were originally compiled by IGruber et al.l ( 1200 lb . 
where the original references can be found. The results of our 
reanalysis of RXTE data are c onsistent with th e earUe r results. 
The INTEGRAL point is from iKIochkov et al.l (l2006l) . The in- 
set at the lower right corner shows the ASM maximum 35 d flux 
(right hand scale). The dotted lines are mean cyclotron line en- 
ergies (Ec) before and after 1991. 



flux (2 - lOkeV), which we take as a measure of the X-ray lu- 
minosity, and the cyclotron line centroid energy Ec. If we take 
the locally measured fluxes (e.g. from our spectral fits) we find 
that they track the maximum ASM flux quite well and they show 
the same correlation with E^. However, we prefer the Main-On 
state maximum flux as the better information on the accretion 
state of the source, since the locally measured flux is subject to 
modulation by variable absorption/shading by the accretion disk. 

For a quantitative analysis of this correlation, which is al- 
ready evident from Fig. 1, we restrict ourselves to our results 
from RXTE data only (see Table 1), because they offer the ad- 
vantage of maximum comparability due to a uniform analysis 
of data from the same instruments. In Fig. 2 we plot the cy- 
clotron line centroid energy Ec as a function of the maximum 
35 d Main-On state RXTE/ASM flux. The dashed line is from 
a linear fit, taking uncertainties of both variables into account 
(see e.g. Numerical Recipies, 15.3), and it defines a slope of 
0.66 + 0.10 (68%) keV/(ASM cts/s). Pearsons linear correlation 
coeflicient is 0.90, the Spearman Rank correlation coefficient is 
0.85 (corresponding to the probabilities to find these correlations 
by chance of < 6 10"^ and < 4 10""*, respectively. We consider 
this correlation as highly significant. 

We also point out here, that the apparent decrease of Ec (and 
the correlated ASM count rate) with time after 1991, seen in 
Fig. 1, is an artifact due to the general variability of the source 
and the random way by which the data happened to be taken, and 
not a secular decrease. Despite strong short-term variability the 
RXTE/ASM light curve of Her X-1 does not provide any indi- 
cation for a long-term (> 5 yrs) change of the mean luminosity. 



4. Variations in cyclotron lines 

Apart from variations with pulse phase, shifts of line positions 
have in phase averaged spectra been found in correlation with 
changing X-ray luminosity in a number of high luminosity tran- 
sient sources. However, these correlations are negative: Ec is re- 
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Fig. 2. The centroid cyclotron line energy of Her X-1 versus 
the maximum flux during the corresponding 35 day Main-On 
as observed by the RXTE/ASM. A linear fit (taking uncertain- 
ties in both variables into account) defines a slope of 0.66 + 
0.10keV/(ASM cts/s). 



duced when Ly increases. iMihara et al.l (Il998l) have interpreted 
this for 4U 0115+63, Cep X-4, and V 0332-H53, as due to a 
change in height of the shock (and emission) region above the 
surface of the neutron star with changing mass accretion rate, 
M . The accretion is in th e super-Eddington regime. In the model 
of iBurnard et al.l (|1991|) , the height of the polar accretion struc- 
ture is tied strongly to M. From this model one expects that an 
increase in accretion rate leads to an increase in the height of the 
scattering region above the neutron star surface, and therfore to 
a decrease in magnetic field strength and hence E^. During the 
2004/2005 outburst of V 0332-H53 a clear anti-correlatio n of the 
line position with X-ray flux was observed (K reykenbo hm et al.l 
120051; iTsvgankov et al.H2006al: iMowlavi et al.ll2006l) . Analysing 
RX TE data of flie Feb-A pr 1 999 outburst of 4U 0115 -H63, 
both iNakaiima et"!!] ( l2006h and iTsvsankov et al.l ( l2006bl) find 
a general anti-correlation between E^ and luminosity. There 
seems to be a steep dependence in a limited luminosity range 
(2...5 • lO^^erg/s), and neighboring regions at lower and higher 
luminosities with near independence of the two quantities. At the 
lowest luminosity the data are even consistent with a reversal of 
the dependence (albeit with low statistical significance). 

For Her X-1, however, so far no clear relation between the 
cyclotron line energy Ec and X-ray flux of any kind had been 
found (see e.g. iMihara et af]|1998l: iGruber et al.ll2001h (but see 
our note in the Discussion). In this Letter we report the first 
detection of a clear positive correlation between E^ and X-ray 
luminosity in Her X-1 (and in any X-ray pulsar, if one re- 
quires a convincing evidence of sufficient significance; see the 
Discussion). 



5. Physics in the sub-Eddington accretion regime 

In the following we suggest a physical explanation of our find- 
ing and attempt a quantitative estimate of the expected effect. 
Since the pioneering paper by iBasko & SunvaevI (Il976l) . it has 
been known that the structure of the accretion column near the 
neutron star surface in X-ray pulsars differs for high-luminosity 
(L > Lc) and low-luminosity (L < L^) regimes, where the criti- 
cal luminosity Lc is given by the local Eddington luminosity Le 
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^ lO^*" erg/s 



(Tj 
CTm 



(1) 



Here M is the NS mass (assumed to be 1.4 Mq in the nu- 
merical estimation), cr„, is the photon-electron scattering cross- 
section in the magnetic field, ctj is the Thomson cross-section, 
and 0c is the half-opening angle of the polar cap magnetic field 
lines. We note that sub-/super-Eddington accretion cannot be 
judged by the observed luminosity alone, decisive is the lo- 
cal accretion rate which depends on the accretion area and the 
photon-electron scattering cross section, which in turn depends 
on the magnetic field (strength and structure). In the radiation- 
dominated super-Eddington regim e the height of the accretion 
column increases with luminosity (iBurnard et al.lll99Tl) . As dis- 
cussed above, this explains the observed anti-correlation be- 
tween the CRSF energy in spectra of bright transient X-ray 
pulsars with luminosity ( | Mihara et al. 1998; Ts vgankov et al.l 
l2006allbl; INakaiima et al.ll2006l) . as the magnetic field value de- 
creases with the height above the NS surface. 

In the sub-Eddington regime of accretion, however, we can 
expect a different behaviour of the CRSF centroid energy when 
the luminosity changes, turning even to the opposite to what is 
observed in high-luminosity X-ray pulsars. For the set of RXTE 
observations of Her X-1 we have a case of sub-Eddington ac- 
cretion in this X-ray pulsar. With the cyclotron line energy E^ 
of about 40 keV, the cross-section of photon-electron scattering 
in the direction along the magnetic field lines (which is of inter- 
est for the radiatio n pressure against th e infalling material) de- 
creases oc (E/Ec)^ jHarding & Lail2006l and references therein). 
At low energies where most of the photons are produced in the 
case of Her X-1, (cr„,lcrj) <K 1, and for an X-ray luminosity of 
a few times lO^*^ erg/s and 6c not much smaller than 0.1 (~ 6°) 
the source is in the sub-Eddington regime. 

The physical picture of magnetic accretion in low-luminosity 
X-ray pulsars was discussed by Nelson et al. (1993). We as- 
sume that the accreting protons lose their kinetic energy in an 
electron-proton atmosphere due to the Coulomb drag and collec- 
tive plasma effects and neglect the possibility for the form ation 
of a collisionless shock (see e.g.lLan ger & RappaportI 19821) . The 
characteristic braking length for protons can be identified with 
their mean free path /, ~ l/riccr, where cr is the effective in- 
teraction cross-section. A change in the observed cyclotron line 
energy can be associated with a change in the height of the emis- 
sion region above the neutron star surface. We write the observed 
cyclotron energy as 



(neB\ ,n 



(2) 



1 /2 

where B is the polar magnetic field strength, g^^ - 
^Jl - Rg/r - 1/(1 + z) is the Schwarzschild metric coefficient 
determining the gravitational redshift z- Rg - IGMjc- is the 
gravitational radius. Identifying the characteristic length Z, with 
the height of the scattering region above the surface of the NS, 
we define r - R + h (with R being the radius of the NS). Then 
the fractional change of the cyclotron line energy is AEc/Ec = 
AB/B + {I /2){Agoo/ goo). Assuming a pure dipole field {B oc r"^) 
we obtain AEc/Ec = -3Ar/R + (l/2){Rg/R)(Ar/R)/gQQ. 

The redshift term can be neglected since for a canonical neu- 
tron star its contribution is only of the order of 10% of the dipole 
term. Making the identification Ar/h - -A«p/«p, we arrive at 



jNelson et al.lll993h 
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Now we need to find how the density in the region of energy 
release and line formation is related to the luminosity. There 
is no simple estimate of this, but we can use the approach by 
iNelson et alj (Il993h to find the structure of the accr etion mound. 
In add ition to the hydrostatic term considered by iNelson et al.l 
(1 1993b . we add the dynamical pressure of infalling protons so 
that the total pressure is P = liiekT - gy + (po^o ~pv^)- Here y is 
the mass column density, T is the temperature, po and vq are the 
density and velocity of the free-falling matter at the beginning 
of the braking region. It is of course not a self-consistent treat- 
ment of the problem, but a step forward to account for dynamical 
pressure of accreting protons inside the braking region. For the 
free fall acceleratio n near the NS surf ace g ^ 2 ■ lO'^^cms"^ 
and y ^ 20 g cm"^ dMiller et al.l 19871) . the hydrostatic term is 
not higher than 4 ■ lO'^dyncm"^. Th e dynamical term ca n be 
expressed using the solution found by INelson et al.l (1 19931) and 
the continuity equation M - pvA (A is the accretion area for one 
pole) in the form po^'o - pv^ = (M/A)vo(l - (1 - r/r,)'''**), where 
T is the Thompson optical depth and t, is the proton stopping 
depth. Clearly, the dynamical term vanishes at the base of the 
accretion column (i.e. at t = t») where the bulk kinetic energy 
of the protons is transformed into their thermal motion. It dom- 
inates, however, at the beginning of the deceleration region (for 
small y and small t, and for typical values like M - 10'^ g s"', 
A-\ /400 of the NS surface area and vq - 10'" crt i s~'). As the 
stopp ing depth is t» ~ 50 for magnetic accretion dNelson et alj 
Il993h . we find for t <s; r» (from where the emission escapes) 
n, - (M I A)(vQl2kT){T I At,). 

As long as /« «: /J we can neglect a variation in A (this is jus- 
tified since the stopping length for r, ~ 50 and n,, ~ 10^"* cm"^ 
is very small), so main parameter which determines the den- 
sity inside the stopping region is M. We also neglect the prin- 
ciple dependence of A on M (due to the change of the disks in- 
ner radius P,„) since this dependence is weak (for dipolar field 
So we finally obtain: 



_ ^ /. AM 



3^^ 
R L 
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We see that in the sub-Eddington accretion regime the frac- 
tional change in cyclotron line energy can be directly propor- 
tional to the fractional change in luminosity. For example, for 
the observed luminosity variations in Her X-1 AL/L ~ 1 (a 
change of a factor of two in the observed flux) the formula above 
gives AEc/Ec ~ 3U/R- The stopping length is /» = T,/{neCrx), 
the density estimation for M - 2 ■ 10'^ g s"', vo = 10'" cm s"', 
kT - lOkeV, and t ~ 1 yields iig ~ lO^cm"-', so h ~ lO^'cm. 
So the expected fractional change in the cyclotron line energy 
is AEc/Ec - 0.03, which is very close to the observed value of 
5%. We note, however, that we used a fiducial value of t, - 50 
appropriate for a weak magnetic field dNelson et alJll993l) . but 
in strong magnetic fields this value can be higher, the electron 
density would decrease and /, increase, accordingly. So we sug- 
gest that a measurement of AE^/Ec can in fact be used to assess 
the stopping length of protons in the regime of sub-Eddington 
magnetic accretion onto neutron stars. 

6. Discussion 

Using standard parameters the calculations above for the sub- 
Eddington accretion regime lead to an estimate which is close to 
the observed effect in Her X-1 (~ 5% in E^ for a factor of two in 
luminosity). The decisive parameter is the local Eddington rate 
at the neutron star Its value depends on the area upon which ac- 
cretion proceeds, so it is expected to vary from pulsar to pulsar. 



In transient pulsars, such as 4U 01 154-63, we may already have 
evidence for a transition from super- to sub-Eddington accre- 
tion ( Tsygankov et al. 2006b; Nakaiima et al. 2006; Terada et aT| 
|2006|) at the decline of the outburst when the luminosity drops 
below ~ 5 ■ 10"*^ erg/s. The sharp jump in the observed cyclotron 
line energy in 4U 01 15-H63 reported bv lTsvgankov et aiTd2006bt) 
could be a sign of an abrupt change in the structure of the accre- 
tion column during such a transition. 

It is conceivable that long-term variations in the accretion 
rate exist also in the persistent X-ray pulsar Her X-1. If so, we 
speculate that the abrupt jump in the cyclotron line energy in 
Her X-1 noticed in the early 1990s (Fig. 1) might be associ- 
ated with a transition of this pulsar from super-Eddington to 
sub-Eddington accretion. We plan to have a close look at the 
historical data with respect to variations in absolute flux (which 
is not an easy task because of the difficulties connected with the 
inter-calibration of different instruments). Taking the binary X- 
ray pulsar A 0535-H262 as an exa mple, the absence of a strong 
luminosity dependence reported by iTerada et alj d2006l) may in- 
dicate that this source stays always in the sub-Eddington regime 
(for example, because the value /» increases with the magnetic 
field strength). If anything , we read a s l ight in crease in E^ with 
luminosity from Fig. 4 of iTerada et al.l d2006l) . Another indica- 
tion of a positive correlation between Ec and l uminosity (and a 
qualita tive interpretation) has been reported by iLaBarbera et al.l 
d2005l) for GX 301-2 (albeit only on the basis of two observa- 
tions and at a significance leve l of about 2 standard deviations). 
While iMihara et all dl998h and lGruber et alj d2001h have shown 
that Her X-1 does not follow the negative correlation known 
from other sources, their data (Table 2 and Fig. 5, respectively) 
do in fact both show a small positive trend (not noticed by the 
two groups of authors at the time). While the statistical signifi- 
cance is low, both sets of data are consistent with our finding of 
5% increase in E^ for a factor of two increase in flux. 

Finally, we note that Her X- 1 is the only highly magnetized 
accreting pulsar for which repeated observations over longer pe- 
riods of time exist. We therefore urge that the source continues 
to be monitored regularly, the source is bright and persistent and 
gives excellent statistics in relatively short observations. 
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